Anomalous Resonance of the Symmetric Single-Impurity Anderson Model in the 

Presence of Pairing Fluctuations 
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We consider the symmetric single-impurity Anderson model in the presence of pairing fluctuations. 
In the isotropic limit, the degrees of freedom of the local impurity are separated into hybridizing and 
non-hybridizing modes. The self-energy for the hybridizing modes can be obtained exactly, leading 
to two subbands centered at ±{7/2. For the non-hybridizing modes, the second order perturbation 
yields a singular resonance of the marginal Fermi liquid form. By multiplicative renomalization, the 
self-energy is derived exactly, showing the resonance is pinned at the Fermi level, while its strength 
is weakened by renormalization. 
PACS numbers: 71.10.Hf, 72.15.Qm, 75.20.Hr 



The unusual normal state of the high-T c cuprate super- 
conductors has been interpreted as a kind of non-Fermi 
liquid (FL) behavior, but so far there have been no fully 
convincing microscopic justifications of such a state in 
quasi-two-dimensional systems with short range interac- 
tions. A phenomenological marginal FL approach was 
suggested by Varma et al. as the "gentlest" departure 
from the conventional FL theory to describe the anoma- 
lous properties. The central hypothesis is that over a 
wide range of momenta the charge and spin polarizabil- 
ities do not contain any intrinsic low energy scales, de- 
pending only on temperature and external frequency: 
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where A^O) is the density of states at the Fermi level, 
and lo c is a high-energy cut-off of the order of 0.5 eV. 
Apart from cuprates, the anomalous behavior observed 
in the normal state of some heavy fermion superconduc- 
tors can also be interpreted in terms of such a marginal 
FL (|] . A number of serious attempts have been made to 
provide a microscopic justification for this phenomeno- 
logical theory, beginning with the three-body resonance 
H , followed by a detailed study of related quantum im- 
purity models |^-^| , and the most recent interpretation is 
based on local fluctuations of the order parameter near 
the quantum critical point using a two-band model for 
cuprates |]. As far as we understand, the issue is still 
open. 

On the other hand, a large body of experimental 
data in NMR, /xSR, specific heat, transport, tunneling, 
neutron scattering, and photoemission measurements in 
most underdoped cuprate superconductors can be inter- 
preted as due to the presence of a pseudo-gap (depleted 
density of states at the Fermi level) above T c [|Io] p"2[ . 
There have been various proposals for the origin of such 
a pseudo-gap, and one of them is due to precursor pairing 
fluctuations p0| . However, a possible connection of such 
pairing fluctuations with marginal FL behavior has not 



yet been explored so far. 

In this Letter, we, following the experience in the heavy 
fermion studies, consider the hybridization of a strongly 
localized electron with an uncorrelated conduction elec- 
tron band, in terms of an impurity model |jl3|| . Encour- 
aged by the recovering of non-FL behavior in the single- 
impurity two-channel Kondo model fl^ , |l5|| , we explore 
a generalized symmetric Anderson model with particle- 
particle mixing on top of the standard particle-hole mix- 
ing. This model can be justified in the presence of pairing 
fluctuations between the local impurity and conduction 
electrons. When the two types of hybridizations have the 
same strengths, the local impurity degrees of freedom are 
separated into hybridizing and non-hybridizing modes, 
and the anomalous polarizability (1) emerges naturally as 
due to the impurity contribution in the non-interacting 
limit. The self-energy for the hybridizing modes can be 
obtained exactly, and the corresponding density of states 
forms two subbands centered at ±{7/2, with U as the 
Hubbard repulsion. Meanwhile, the self-energy for the 
non-hybridizing modes is calculated in the second order 
perturbation, giving rise to a singular resonance at the 
Fermi level, of the marginal FL type. Using the mul- 
tiplicative renormalization method, the exact self-energy 
can be derived self-consistently in the weak coupling limit 
and extended then to strong coupling. We find that the 
non-FL resonance is pinned at the Fermi level and it is 
described by the X-ray edge type of singularity. 

The symmetric single impurity Anderson model in the 
presence of pairing fluctuations is given by: 

iJ = Ve,-Ct C s + -V(C*t d a +dlC r ) 
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where the symmetric condition e c i = —{7/2 has been as- 
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sumed, and the chemical potential of conduction elec- 
trons is set to zero. It is worthwhile to note that the 
V a term is similar to the transverse interaction term for 
the second channel in the effective two-channel Kondo 
Hamiltonian obtained in the abelian bosonization ap- 
proach |Q. We will focus on the case of V a = V, cor- 
responding to the isotropic case of the model. As jus- 
tification of this model we would mention that Varma 
has emphasized the importance of the potential scatter- 
ing at the impurity Vq J2k,k',a c k,a c k' ,<rd\d a @7 an d the 
renormalization group analysis showed that the unitar- 
ity limit for this potential scattering is being approached, 
while the opposite-spin scattering is renormalized to zero 
. A mean field decoupling of this quartic term allowing 
"anomalous average" {d a c a ) due to pairing fluctuations 
would lead to the particle-particle mixing term in Eq. 
(2). Such a pairing is of triplet character, and whether it 
is related to the "odd frequency" pairing discussed earlier 
fl7f is not considered here. Nor we are going to discuss 
the self-consistency issue of the pairing per se. We will 
rather focus on the physical consequences of the above 
well-defined model. 

In order to reveal the characteristics of this gener- 
alized model, it is very useful to introduce the Majo- 
rana fermion representation fl5| : = [di — £d2)/v2, 
di = (— e?3 — i(io)/v / 2. The model Hamiltonian can then 
be rewritten as H = Hq + Hi, where 



u /L^i k k,a k,i 



k,a 
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Now it is clear that the pairing fluctuation terms have 
changed the structure of the model: one half degrees 
of freedom of the local impurity {d\,da) hybridize with 
the conduction electrons, while another half degrees of 
freedom (do,^) decouple. The former will be referred 
to as hybridizing modes, whereas the latter as non- 
hybridizing modes. These two types of modes are cor- 
related through the Hubbard interaction. We emphasize 
that this change is non-perturbative in nature and makes 
the non-interacting limit nontrivial. 

In the U — limit, Hq is solved exactly. The 
Green's functions for the hybridizing and non-hybridizing 
modes are calculated as G[°'(u n ) = - — -r-^ — ? — % and 

G nlM = ^, where u n = (2n + l)n/0 with = 1/T, 
A = irp(0)V 2 is the hybridization width, p(0) is the 
density of states of conduction electrons at the Fermi 
level, and A will be the high energy cutoff in the weak 
coupling theory. The Green's function for the non- 
hybridizing modes describes a fermionic zero mode and 
— — sgn(r)/2. The total impurity density of 



states is thus obtained as Aj(u>) 



1 A 

TT Ul 2 +A 2 
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Apart from a Lorentz distribution, a delta function peak 
due to the non-hybridizing modes appears at the Fermi 
level. As a result, the impurity charge and spin density 
fluctuations will be affected. Since the charge and spin 
density operators are defined as 



1 i 

n d = ^(<4 d T + d \ d l ~ 1 ) = -9(^2 + d d 3 ), 
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(4) 



the charge and spin dynamical susceptibilities are both 
given by 

= 4 0) K)^(^ - (5) 

where VL n = 2nn/0. After summation over frequency, 
the retarded imaginary part is obtained: 



ImxW(n,T)=-|iV(0)tanh^ 



(6) 



where N(fY) — 7 02+ a 2 is the density of states of the 
hybridizing modes. This is exactly what has been as- 
sumed in the marginal FL phenomenology jjj, of the 
same form as Eq.(l)! It becomes transparent that the 
anomalous impurity charge and spin polarizability is due 
to the delta resonance at the Fermi level, which is induced 
by the pairing fluctuations under the condition V a = V. 
An analogous situation previously appeared in the stud- 
ies of the two-channel Kondo model in the isotropic case 

Now consider the effects of interaction on the impurity 
density of states. The perturbation theory in terms of 
Majorana fermions has been fully developed in the earlier 
publications Gq]. In the second order perturbation, the 
self-energy of the hybridizing modes is given by: S^(t — 
u 2 M0), 



4T"h ( T 



t'). It is worthwhile to note that this 
self-energy is exact, which can also be obtained by means 
of retarded double-time Green's functions with the help 
of [d^da, H] — 0. So the exact Green's function is 

G^K) = [Gi 0) K)] 1 ^G< 0) K), (7) 
and the density of states ImG^(w) is thus 



1 

2^ 
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(u - U/2) 2 + A 2 (w + 172) 2 + A 2 



(8) 



This spectral density is similar to that of the conven- 
tional symmetric single impurity Anderson model in the 
Hartree-Fock approximation |13|| . 

Meanwhile, the self-energy of the non-hybridizing modes 
in the second order perturbation is expressed by: 



tr 2 , 
" p 2 ■ 



Gf ( Wl )G^(o; 2 )G^K (9) 
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Using the spectral representation, the hybridizing Green's 
functions are expressed in terms of their density of states 
N(e), and after the frequency summation we find 



N(e 1 )N(e 2 ) 



- ei 



tanh (£) 



tanh ( — ^ ) 
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Carrying out the analytical continuation iuj n — ► cj and 
approximating iV(e) by its value at the Fermi level, the 
imaginary part of the retarded self-energy is derived as 



ImE$(w,T) 



7rA 



' coth 



-1 
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and the corresponding real part is also obtained as 

ReE^T) « (J^) 2 w ln(2^) . This self-energy 

gives rise to a quasiparticle weight depending logarithmi- 
cally on frequency or temperature, namely, 



,(2) 
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(12) 



In the corresponding Green's function, the pole u> = 
is preserved and the resonance is characterized by a 
marginal FL origin. To some extent, the non-hybridizing 
modes may represent the low energy excitations of the 
model. 

Since there is a logarithmic correction to the self- 
energy for the non-hybridizing modes, we have to ex- 
amine the higher order expansions carefully. Moreover, 
the logarithmic singularity in the reducible interaction 
vertex diagrams can be proved to cancel each other, so 
one has to consider the irreducible interaction vertex 
ri,2,3,o(0, u> n ; 0, oj n ) = T(u n ), which can not be severed 
into separate diagrams by cutting a pair of lines, one of 
which is the propagator of G n h and the other is Gh- The 
lowest order correction to T(uj n ) is given by 



U 3 



GfK^rM Gnhfan - Wl - W 2 ) 
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(0), 



It is interesting to note that this lowest order correction 
can be related to E^?(u; n ) by the relation r( 3 )(w„) = 

— JJ d a?' 1 ^"' 1 ■ Therefore, the correction to the irreducible 
vertex is also logarithmically divergent: 



rW( Wl T)«-cr(^ 



In 
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(13) 



In fact, from a general consideration of the model, a Ward 
identity can be established relating the irreducible 
vertex to the non-hybridizing self-energy: 



r(« n ) = u 



1 



d(iuj n ) 



(14) 



which shows that to a large extent the singularity in 
r((x>n) and Ti n h{to n ) cancel each other. This implies that 
the renormalization of the self-energy in the higher order 
expansions must be treated on an equal footing with that 
of the irreducible vertex. 

In the treatment of logarithmic problems, the multi- 
plicative renormalization is an effective method pOj . In 
such a treatment, the temperature is fixed at zero for 
simplicity and only the frequency variables are retained. 
Since the interaction is cut off at A in our weak coupling 
theory, the Green's functions and vertices will depend 
only on the ratio ^, and the dimensionless coupling con- 



stant U 



must be defined as the proper invariant 



(11) coupling constant. When we express the irreducible ver- 



tex by = UT, it follows from the structure of the 
Dyson equation that the renormalization equations will 
have the forms 

G h (w,A',U') = z 1 G h (w,A,U), 
G nh {u/A',U') = z 2 G nh (uj/A,U), 
f(u;/A',U') = z- 1 f(u J /A,U), 
U' = z^z^zaU. (15) 

Since its second order perturbational result is exact, the 
Green's function for the hybridizing modes should remain 
unrenormalized and therefore z\ — 1. Using the low order 
perturbation results for z n h and T, we deduce 



A' 



A' 



= 1 + {/ In I — ) , ,z 3 = 1 — U'hx I — I , (16) 

and U' = U. The invariant coupling constant is not 
renormalized, and this conclusion is actually valid be- 
yond the multiplicative renormalization approximation, 
because it is ensured by the Ward identity! 

Meanwhile, the quasiparticle weight z n h(uj) is also a 
proper quantity satisfying the criterion of multiplicative 
renormalization, and has good transformation properties. 
When the second order perturbation result of z n h is taken 
into account, the Lie equation up to the second order can 
be derived as follows: 



d In Znh _ jy 2 1 
d In bj uj 



(17) 



By integrating over ui and determining the constant of 
integration from fitting to the second order perturbation 

u 2 

expression, we get z n h{uj) = (^) . The Green's func- 
tion for the non-hybridizing modes is thus obtained 



(18) 



To calculate the spectral density lmG n h we note that 
the branch cut of G n h extends only on one side of the 
branching point. Thus, ImG n /j is a step function of the 
energy. The phase of the factor uj a is then completely 
determined: on the branch cut side it involves a factor 
(—1)". In the weak coupling limit U « 1, we thus obtain 
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ImGnhM =*U 2 ^ (£) U ■ (19) 

The singularity of the resonance at to = is weakened 
by renormalization, and follows a power law with an ex- 
ponent proportional to the squared interaction strength. 
The non-FL quasiparticlc weight vanishes at the Fermi 
level, thus there are strong similarities between the 
present singular resonance and the X-ray emission spec- 
tra pH , an d the marginal FL resonance is actually de- 
scribed by a power law behavior in the low energy limit. 
Paralleling the asymptotically exact calculation for the 
X-ray edge problem [p2[ , we believe that it is possible 
to develop an exact treatment of this resonance for arbi- 
trary coupling parameter U. This way, the above weak 
coupling results may be extended to strong coupling limit 
after replacing U by S/tt, where S = arctan([//A) is an 
s-wave phase shift at the Fermi level. 

In addition, when the large U limit is considered, 
a Schricffcr-Wolf transformation can be applied to the 
model Hamiltonian, generating an s-d type of model: 

4V 2 

H ' = E HC{C^ + —l*v(0) + t"(0)]S* (20) 

k,cr 

where <r y (0) = -i(Cl^C 0tl - C^ ; C ,t)/2, and t»(0) = 

— z(Cq |Cq i — Co ) jCo ) f)/2, are the y-components of the 
spin and isospin densities of the conduction electrons at 
the impurity site. S% — —id^di is the corresponding 
impurity spin operator. This resulting model is remi- 
niscent of the so-called compactified two-channel Kondo 
model However, here is only one component in 

the exchange interactions, which become a genuinely 
marginal operator, being compatible with the renormal- 
ization analysis of the Hubbard interaction. 

The total density of states of the local impurity is 
a sum of the hybridizing and non-hybridizing modes: 
ImG/j-t- ImG n h- Within the present exact theory, a very 
interesting picture emerges in a nature way: The hy- 
bridizing modes form two subbands with Lorentz form 
centered at ±?7/2, in some sense similar to the lower and 
upper Hubbard bands, while the non- hybridizing modes 
form a singular non-FL resonance at the Fermi level. The 
low-energy excitations and high-energy excitations are 
thus separated, and only the former is strongly renor- 
malized. More interestingly, a very similar picture with 
three peaks in the density of states has been suggested 
by Kotliar, Georges, and coworkers ]23| in their studies 
of the symmetric d = oo one-band Hubbard model close 
to the metal-insulator transition, where the central peak 
is a FL quasiparticle resonance. It is very likely that 
the present model including the pairing fluctuations be- 
tween the local impurity and conduction electrons could 
shed some light on the nature of some strongly correlated 
electron systems. 

In conclusion, we have considered the effect of pairing 
fluctuations on the symmetric single impurity Anderson 



model. The main features of this model are very simi- 
lar to those anticipated for the "three body" resonance 
model, considered earlier ||. However, the difficulty due 
to generating a new low-energy scale, inherent in the pre- 
vious model, does not appear here. This model study 
suggests a possible microscopic origin of the marginal FL 
behavior. After completing the present calculation, we 
saw a new preprint fl24| where an attempt was made to 
make connection between the marginal FL behavior and 
a lattice of three-body bound states. However, the issue 
considered there is different from ours. 
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